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Abstract. - The age and growth of the moray eel ( Muraena helena) from the north coast of Tunisia were deter¬ 
mined from otolith growth ring studies. A total of 221 specimens were used. Total length (TL) and weight (TW) 
ranged from 31 to 100 cm and from 45.4 to 2311.9 g, respectively. The transverse section of the otoliths observed 
under transmitted light proved to be a useful method for ageing M. helena. Marginal increment analysis was 
used to validate the periodicity of increment formation. The age of the specimens ranged from 2 to 10 years in 
males, and from 2 to 9 years in females. The von Bertalanffy growth parameters were Loo = 198.1cm, k = 0.063 
year 1 2 and to = -0.797 years for the combined sexes. Length-weight relationship was determined as follows: 
W = 0.0013TL 3084 (r = 0.948). 


©SFI 

Received: 18 Aug. 2015 
Accepted: 2 Feb. 2016 
Editor: G. Duhamel 


Key words 

Muraenidae 
Muraena helena 
Mediterranean Sea 
North coast of Tunisia 
Age 
Growth 


Resume. - Age et croissance de Muraena helena (Muraenidae) des cotes nord de la Tunisie. 

L’age et la croissance de Muraena helena ont ete etudies a partir de 221 specimens collectes sur les cotes 
nord de la Tunisie, dont 104 femelles et 117 males. Leur longueur totale variait de 31 a 100 cm pour des masses 
totales respectives de 45,4 a 2311,9 g. L’age a ete estime en comptant les anneaux sur des coupes transversales 
des otolithes (sagittae). L’analyse de Tallongement marginal a mis en evidence un ralentissement en hiver, cor- 
respondant a la formation de l’anneau d’arret de croissance. La determination de Page nous a permis de distin- 
guer neuf classes d’age pour les males et huit pour les femelles. L’equation de croissance de Von Bertalanffy, 
etablie pour les sexes combines, presente les parametres suivants : Loo = 198,1 cm, k = 0,0632 an' 1 , to = -0,797 an. 
Pour les sexes regroupes, la croissance en masse de M. helena est isometrique par rapport a la taille du poisson 
(b = 3,084). Chez les males la masse et la taille augmentent suivant une allometrie majorante (b = 3,229) alors 
qu’elle est isometrique pour les femelles (b = 2,894). 


The family Muraenidae includes 197 species (Eschmey- 
er and Fong, 2015). In the north-eastern Atlantic Ocean and 
the Mediterranean Sea, this family is represented by three 
demersal predator species: Muraena helena Linnaeus, 1758, 
Gymnothorax unicolor (Delaroche, 1809) and Enchelycore 
anatina (Lowe, 1838). The latter was first recorded in the 
eastern Mediterranean by Ben-Tuvia and Golani (1984), and 
now extends west towards Malta (Deidun et al., 2015). The 
Mediterranean moray eel (M. helena) is widely distributed, 
from the British Isles to Senegal (Bauchot, 1986), including 
the Canary Islands, Madeira, and the Azores (Randall and 
Golani, 1995), the two basins of the Mediterranean Sea, east¬ 
ern and western, and the Adriatic Sea (Fisher et al ., 1987). It 
has also been reported off the north coast of Tunisia by sev¬ 
eral authors (Gruvel, 1926; Bourgeois and Farina, 1961; Ben 
Mustapha, 1966; Azouz, 1971,1974; Bradai etal., 2004). 

Because of its aggressive behaviour and its habit of hid¬ 
ing in cracks and crevices, the capture of M. helena remains 
difficult, and studies about the biology of this species are 
scarce. Only a few studies have been published about its age 
and growth (Matic-Skoko et al., 2011 for the Adriatic Sea), 


food preferences (Matic-Skoko et al., 2014) or morphologi¬ 
cal and molecular analyses (Jimenez et al., 2007). However, 
its large distribution and abundance in certain parts of the 
Mediterranean make it an important species, for which it 
is important to obtain reliable information on its age and 
growth, especially for population dynamics studies. 

In teleost fish, age estimation is generally based on scle- 
rochronology, i.e. the study of growth marks deposited in 
hard tissues: bones, scales, and otoliths (Panfili etal., 2002). 
Because M. helena is almost scaleless, we used otoliths (sag¬ 
ittae) to study the relation between fish size and age, and to 
assess growth parameters in moray eels captured along the 
north coast of Tunisia. 


MATERIALS AND METHODS 

Landings from three places (i.e., Bizerte, Kelibia and El 
Haouaria) along the northern coast of Tunisia (Fig. 1) have 
been studied between December 2011 and December 2013. 
A total of 221 M. helena specimens were collected from lon- 
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Table I. - Number of Muraena helena individuals sampled per month between December 2011 and December 2013 from the north coast of 
Tunisia (cumulative data). 


Month 

December 

January 

February 

March 

April 

May 

June 

July 

August 

September 

October 

November 

Number of 
specimens 

23 

31 

14 

28 

12 

10 

23 

20 

15 

10 

14 

21 



Figure 1. - Map of the study area, indicating the three sampling 
locations Q¥) (North coast of Tunisia, central Mediterranean). 


gline and bottom trawlers, caught at depths ranging from 35 
to 200 m. A mean of approximately 18 individuals (range 10 
to 31) were examined per month (Tab. I). Length and weight 
ranged from 31 to 100 cm, and 45.4 to 2311.9 g, respectively. 
For each sample, sagittae were extracted, dried, and stored 
for future study. The removal of the otoliths was difficult 
since M. helena is characterized by a highly ossified skull. 
Sex was recorded after a macroscopic and/or a microscopic 
observation of the gonads. All specimens were sexed and the 
sex ratio (Sr), percentages of male (%M) and females (%F) 
were determined according to: 

Sr = (number of males/number of females) 

%M = percentage of masculinity = (number of males/number 
of males and females) x 100 

%F = percentage of femininity = (number of females/number 
of males and females) x 100 

The sex ratio was analysed using the Chi-square test 
(X 2 )- 

The treatment of otoliths was performed at the platform 
of otolithometry, in the laboratory of ichthyology of the 
Museum national d’Histoire naturelle, Paris. 


Sections of sagittae were used for age determination. The 
sagittae were embedded in a polyester resin then sectioned 
transversely through the nucleus and polished to obtain 
transverse sections without scratches. 

Sections were mounted on glass microscope slides and 
examined for banding patterns under a binocular microscope 
(x40) with transmitted light to detect annuli. A drop of rose¬ 
mary essential oil was placed on the sections to enhance the 
contrast of rings. All otolith sections were read independ¬ 
ently by two researchers. 

Marginal increment (MI) analysis was carried out to 
validate the periodicity of annulus formation in otoliths. MI 
was calculated as the distance between the outer edge of the 
last annulus to the edge of the otolith (Hood et al ., 1988; 
O’Sullivan et al ., 2003). 

For simulating growth in M. helena , we used the von 
Bertalanffy growth function (VBGF), defined by : Lt = Loo (1 
- e _k(t_t °)), where Lt is the total length at age t; Loo (cm) is the 
asymptotic length; K (year -1 ) is the rate at which the asymp¬ 
totic length is approached; t (year) is the age of the fish; and 
to (year) is the hypothetical age at zero length. The von Ber¬ 
talanffy growth parameters (Loo (cm), k and to (years)) were 
calculated for the whole sample using the method of Ben 
Salem and Daget (1991). 

For each sample, total length (TL, cm) and total weight 
(TW, g) were recorded. TL was measured to the near¬ 
est 0.1 cm, and TW to the nearest 0.01 g. The relationship 
between weight and total length is not linear and is best 
described by a power function of the form W = a TL b , where 
W = weight in grams, a = coefficient factor, TL = total length 
in centimeters, and b = allometry coefficient. In the most 
common pattern, the mass varies as the cube of the length 
(Ricker, 1975; Busacker et al., 1990). An allometric coeffi¬ 
cient b value shows an allometric growth when it is larger 
or smaller than 3, and an isometric growth when it is equal 


Table II. - Size (Lt) and weight (W) for the sampled specimens of Muraena helena. n: sample size; m: mean; min: minimum value; max: 
maximum value; SD: standard deviation. 



Male (n = 117) 

Female (n = 104) 

Combined data (n = 221) 


m 

min 

max 

SD 

m 

min 

max 

SD 

m 

min 

max 

SD 

Lt (mm) 

720.4 

330.0 

1000.0 

128.5 

676.5 

310.0 

925 

119.8 

699.7 

310.0 

1000.0 

126.1 

W(g) 

785.6 

60.9 

2311.9 

471.6 

599.6 

45.4 

1780 

302.4 

698.1 

45.4 

2311.9 

410.8 


Table III. - Numerical proportions of males, females and overall sex-ratio of Muraena helena from the north coast of Tunisia, n: sample 
size; %M: percentage of masculinity; %F: percentage of femininity; Sr: sex ratio; x 2 - chi-square test. 


Period 

n 

n male 

n female 

%F 

%M 

Sr 

X 2 

P 

Significance 

Dec. 2011-Dec. 2013 

221 

117 

104 

47.06 

52.94 

1.125 

0.764 

0.382 

not significant 
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to 3. In order to confirm whether the b value obtained in the 
power regressions was significantly different from the iso¬ 
metric value (b = 3) the t-test was used. 


RESULTS 

The descriptive statistics of total length and weight for 
the 221 specimens are presented in table II for both sexes and 
for combined data. Among the 221 specimens of M. helena 
sampled, 104 were females and 117 were males. The male 
to female ratio was 1.12:1.00, which does not differ statisti¬ 
cally from the expected 1:1 (Tab. III). The length frequency 
distribution of the total sample indicates that the bulk of the 
catch were between 60 and 80 cm TL (Fig. 2). 



TL (cm) 


Figure 2. - Length-frequency distribution of females, males and 
combined sexes for the Muraena helena sampled from the north 
coast of Tunisia. 



20[jm 

Figure 3. - Transverse section of an otolith (sagitta), observed under 
a binocular microscope with transmitted light. Mureana helena 
from the north coast of Tunisia (TL = 62.3 cm, W = 442.73 g). 
Black spots indicate the opaque bands. 


Ageing and validation 

Under transmitted light, otolith sections showed a dark 
nucleus surrounded by alternating translucent (annuli) and 
opaque bands (Fig. 3). The formation of growth increments 
was validated through marginal increment analysis. The 
age of each fish was attributed by counting the number of 
annuli. The analysis of seasonal variations in the mean mar¬ 
ginal increment showed that the formation of growth incre¬ 
ments clearly followed a seasonal pattern (Fig. 4), and indi¬ 
cates that the mean MI presents a minimum value in winter 
(0.059), resulting in the appearance of a growth arrest ring. 
Assuming that one annulus was deposited regularly each 
year, the counts revealed 9 age classes, ranging between 2 
and 10 years. The oldest males and females of M. helena 
were estimated to be 10 and 9 years old, respectively. Age 
class 7 predominated (26.2% of individuals) in both males 
and females within the total catch sample. Interestingly, age 
classes 6 and 7 represented 46.6% of individuals in the total 
sample, a frequency very close to that found by Matic-Skoko 
et al. (2011) in the Adriatic Sea. 

We grouped individuals according to the number of 
translucent bands they presented, which allowed us to assign 
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Figure 4. - Seasonal mean marginal increment observed from sagit- 
tae sections of ages 2 to 10 Muraena helena specimens. 



Figure 5. - Relationship between age and total length of Muraena 
helena (n = 221) sampled from north coast of Tunisia. 
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Table IV. - Age-length key for all specimens of Muraena helena from the north coast of Tunisia, m: mean; SD: standard deviation; TL: total 
length; W: weight. 



Age groups (years) 

Length interval 
(cm) 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Total 

30.1-32.5 

2 









2 

32.6-35.0 

2 









2 

35.1-37.5 

1 

1 








2 

37.6-40.0 










0 

40.1-42.5 










0 

42.6-45.0 


2 








2 

45.1-47.5 



1 







1 

47.6-50.0 



3 







3 

50.1-52.5 



4 

2 






6 

52.6-55.0 



3 

1 






4 

55.1-57.5 



4 

1 






5 

57.6-60.0 



3 

6 






9 

60.1-62.5 



2 

10 

3 





15 

62.6-65.0 




15 

13 





28 

65.1-67.5 




5 

14 





19 

67.6-70.0 




1 

10 

2 




13 

70.1-72.5 





2 

14 




16 

72.6-75.0 





2 

15 

3 



20 

75.1-77.5 






13 

1 



14 

77.6-80.0 





1 

6 

3 



10 

80.1-82.5 






3 

12 



15 

82.6-85.0 






5 

12 



17 

85.1-87.5 







4 



4 

87.6-90.0 







1 

2 


3 

90.1-92.5 







1 

2 


3 

92.6-95.0 








2 


2 

95.1-97.5 








2 

1 

3 

97.6-100.0 








1 

2 

3 

Total 

5 

3 

20 

41 

45 

58 

37 

9 

3 

221 

% 

2.26 

1.35 

9.04 

18.55 

20.36 

26.24 

16.74 

4.07 

1.35 

100 

m TL (cm) 

33.28 

41.06 

54.12 

61.95 

66.68 

75.40 

82.10 

93.58 

99.00 


SDTL 

21.89 

35.23 

43.48 

38.28 

33.20 

37.95 

37.59 

34.48 

17.32 


m W (g) 

64.56 

109.20 

292.63 

430.73 

541.01 

834.29 

967.78 

1688.94 

2122.95 


SDW 

19.62 

38.64 

88.91 

76.17 

152.73 

171.25 

284.80 

505.67 

5.80 



Table V. - Growth parameters for males and females of Muraena 
helena from the north coast of Tunisia. L^: asymptotic length; k: 
growth coefficient; to: theoretical age at zero length. 



Loo (cm) 

K (yr 1 ) 

to (yr) 

Males (n = 117) 
Females (n = 104) 

199.3 

194.3 

0.064 

0.059 

-0.654 

-1.3 


an average length for each age group. Table IV summariz¬ 
es length and age data for the M. helena sampled from the 
north coast of Tunisia, while figure 5 shows the relationship 
between TL and age. Despite some overlapping of indi¬ 
viduals of different lengths corresponding to the same age, 


and vice versa, total length and age were well correlated 
(R 2 = 0.9). 

The age frequency distribution of M. helena (Fig. 6) 
shows that females were better represented in the age classes 
2, 3, 4 and 5, with respective percentages of 2.9, 1.9, 10.6 
and 21.1%. In older age groups, especially in the age group 9 
and 10, males were more numerous. 

Growth model - Length-weight relationship 

The VBGF for the entire population sampled had the 
following parameters: Loo = 198.1cm, K = 0.063 yr" 1 and 
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Figure 6. - Age distribution of males (M), females (F) and combined 
sample (T) of Muraena helena from the north coast of Tunisia. 



Figure 7. - Length-weight relationship in Muraena helena sampled 
between December 2011 and December 2013 (north coast of Tuni¬ 
sia). 


t 0 = -0.797 yr. The estimated parameters for males and 
females are provided in table V. Growth was similar in both 
sexes, although females grow at slightly slower rate than 
males. 

The calculated length-weight equation for the whole sam¬ 
ple was W = 0.0013 x TL 3 - 084 (Fig. 7). The estimated param¬ 
eters of the length-weight relationship for the total sample 
of males and females are shown in table VI. Males showed 
positive allometric growth with a b value (3.229) significant¬ 
ly different from 3. However, the values of b for females and 
combined sex were 2.894 and 3.084, respectively, indicating 
an isometric growth. A correlation coefficient close to 1, for 


the total sample, male and female, showed that the growth of 
the fish mass was highly correlated with its size (Tab. VI). 

DISCUSSION 

Age data of fish populations is highly important, as the 
majority of biological estimations in fisheries are based on 
this parameter (Campana, 2001). Otolith reading is a long 
established method in fisheries for ageing teleosts (e.g., 
Hederstrom, 1759; Reibish, 1899; Cunningham 1905; Wal¬ 
lace 1907), with interest increasing considerably since the 
discovery of daily growth marks by Pannella (1971). How¬ 
ever, Muraena helena is still a poorly studied species, for 
which only two published data sets are available for its age 
analysis and growth: one for the Canary Islands (Jimenez 
et al ., 2007), the other for the Adriatic Sea (Matic-Skoko et 
al., 2011); both of which used the observation of otoliths to 
determine age. The sagittae of M. helena clearly exhibited 
strong sclerochronological information, with cyclical growth 
marks (zones and annuli) easily identifiable and readable on 
transverse sections, confirming the results of Matic-Skoko et 
al. (2011). 

It is impossible to estimate accurately individual age 
without precisely identifying the growth marks and the tim¬ 
ing of their formation (Panfili and Ximenes, 1994). In an 
attempt to validate the annual deposition of bands in the 
otolith (Campana, 2001), marginal increment analysis was 
carried out. We determined that the opaque zones are com¬ 
pletely formed during the summer, while translucent zones 
are deposited during the autumn/winter period. The mean 
marginal increment was lowest during the winter; a season 
characterized by a fall of temperatures and a lack of prey, 
two factors impacting severely fish growth. Moreover, since 
M. helena rarely leaves its territory, except for short distanc¬ 
es in search of prey (Bohlke et al ., 1989), it is particularly 
sensitive to changes in local environmental conditions. 

According to our results, the age of M. helena from the 
commercial catch, along the northern coast of Tunisia, var¬ 
ied from two to 10 years. Matic-Skoko et al. (2011) indicat¬ 
ed that the age of this species varied from two to 12 years in 
the Adriatic Sea, while Jimenez et al. (2007) reported that 
the age of M. helena in the Canary Islands varied from three 
to 15 years. It is noteworthy that during our survey no speci¬ 
mens under two years old were captured, maybe because of 
fishing selectivity or because they 
are not yet sedentarized. The max¬ 
imum age observed on the north 
coast of Tunisia was 10 years for 
a male of 100 cm total length. It is 
possible that individuals over 10 
years exist in Tunisia as the maxi¬ 
mal size recorded in this study 


Table VI. - Length-weight relationships (W = aL b ) for males, females and total sample of 
Muraena helena from the north coast of Tunisia during December 2011-December 2013 period. 
R: coefficient of correlation; N: number of specimens. 



Equations 

b 

R 

N 

Lai (to .05 - 1-98) 

Allometry 

Males 

Wt = 0.0007 TL 3 229 

3.229 

0.960 

117 

2.600 

Positive allometry 

Females 

Wt = 0.0027TL 2894 

2.894 

0.931 

104 

-0.844 

Isometry 

Total sample 

Wt = 0.0013TL 3084 

3.084 

0.948 

221 

1.187 

Isometry 
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Table VII. - Parameters of the length-weight relationship of Muraena helena reported by authors in 
three localities, b: allometric coefficient; R 2 : coefficient of determination. 



b 

R 2 

Allometry 

Localities 

Authors 

Sexes combined 

3.314 

0.887 

Positive 

Canary Islands 

Jimenez et al., 2007 

Sexes combined 

2.776 

0.817 

Negative 

Southern Adriatic Sea 

Matic-Skoko et al., 2011 

Sexes combined 

3.084 

0.898 

Isometric 

North coasts of Tunisia 

Present study 


Table VIII. - Biogeographic comparison of Von Bertalanffy growth function parameters of Muraena helena. n: number of specimens; Loo: 
asymptotic length; K: growth coefficient; to: theoretical age at zero length. 


Sex 

Age 

(year) 

n 

Age 

determination 

Ltmin L t max 

(cm) 

Loo (cm) 

K (an 1 ) 

to (an) 

Area 

Reference 

Sexes 

combined 

2-12 

193 

Otoliths 

27.5-121.0 

162.7 ±0.489 

0.089 ±0.003 

-0.660 ±0.219 

Southern 

Adriatic 

Sea 

Matic- 
Skoko et al. 
(2011) 

Sexes 

combined 

3- 15 

179 

Otoliths 

49.0-123.0 

170.0 

0.078 

0.355 

Canary 

Islands 

Jimenez et 
al. (2007) 

Sexes 

combined 

2-10 

221 

Otoliths 

31.0-100.0 

198.1 

0.0632 

-0.797 

North 
coast of 
Tunisia 

Present 

study 


(100 cm) is substantially lower than in the two above cited 
studies (121 and 134 cm, respectively). The absence of large 
specimens might be due to fishing pressure, as moray eels 
are commonly found on the Tunisian fish markets, or to other 
reasons (e.g., selectivity of fishing gear, depth). 

Mediterranean moray eels from the north coast of Tuni¬ 
sia showed positive allometric growth for males and iso¬ 
metric growth for females. Length-weight relationships can 
provide useful information about the increase in weight of a 
population, and this parameter could also be important for 
comparative studies between populations. Matic-Skoko et 
al. (2011) reported a negative allometric growth (b = 2.776) 
in the combined sex of M. helena collected from the south¬ 
ern Adriatic Sea. Nevertheless, the b value determined by 
Jimenez et al. (2007) for M. helena from the Canary Islands 
was 3.314, indicating positive allometric growth. The length- 
weight relationship parameters for M. helena from various 
areas are given in table VII. 

Comparison of growth parameter estimates of M. helena 
from different geographical areas showed slight discrepan¬ 
cies (Tab. VIII), but global VB models were similar. The 
biogeographical differences between these regions could 
possibly explain these growth rate variations, as it has 
already been observed in the latitudinal populations of sev¬ 
eral species (e.g., Parrish et al ., 1985). Along the north coast 
of Tunisia M. helena has a lower growth coefficient (k) and 
a higher estimated asymptotic length than those reported for 
the Adriatic Sea or the Canary Islands. M. helena is a slow 
growing and long living Anguilliformes known to attain a 
large size, up to 2 m, which is in agreement with the asymp¬ 
totic length of our VB model (198.1 cm). 


CONCLUSION 

Our study confirms that otoliths are valid and useful age¬ 
ing structures for moray eels. These new data on age and 
growth of M. helena from Tunisian waters have enhanced 
our knowledge of this species in the Mediterranean. How¬ 
ever, the life cycle of Mediterranean moray eel remains 
poorly understood, with further biological studies needed to 
enable better resource management in Mediterranean coastal 
zones, where it has been exploited since at least the Roman 
period. In particular, several points require further research: 
the recruitment and growth of juveniles (< 2 + ), and the exact 
mechanism for geographical differences in growth rates of 
Mediterranean moray eels; since, according to Jimenez et 
al. (2007), M. helena shows little genetic divergence among 
specimens from different localities. 
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